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Purpose: Allergic asthma (AA) and rheumatoid arthritis (RA) are immune tolerance-related diseases, and immune tolerance is known to be influ- 
enced by costimulatory molecules. In this study, we sought to identify common genetic susceptibility in AA and RA. Methods: Two hundred cases 
of AA, 184 cases of RA, and 182 healthy controls were recruited at the Seoul National University Hospital, Seoul, Korea. Eight single nucleotide 
polymorphisms (SNPs) in five genes coding costimulatory molecules, namely, -3180T, +49A>G, and 6230G>A in CTLA4, IVS3+17T>C in CD28, 
-3479T>G and 1179V in CD8B, -1C>T in CD40, and -3458A>G in CD40LG were scored, and genetic interactions were evaluated by multifactor di- 
mensionality reduction (MDR) analysis. Results: MDR analysis revealed a significant gene-gene interaction between -3479T > G CD86 and -3458A> 
G CD40LG for AA. Subjects with the T/T genotype of -3479T>G CD86 and the A/A genotype of -3458A>G CD40LG were found to be significantly 
more likely to develop AA than those with the T/T genotype of -3479T>G CD86 and A/- genotype of -3458A>G CD40LG (adjusted OR, 6.09; 95% 
CI, 2.89-12.98; logistic regression analysis controlled by age). Similarly those subjects showed a significant risk of developing RA (adjusted OR, 
39.35; 95% CI, 15.01-107.00, logistic regression analysis controlled by age). Conclusions: Our findings suggest that a genetic interaction be- 
tween CD86 and CD40LG favors the development of both AA and RA. 
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INTRODUCTION 

Immune tolerance refers to a specific non-reactivity to a given 
antigen which under other circumstances would possibly induce 
an immune reaction. The mechanisms underlying immune tol- 
erance are complex and are influenced by several factors, in- 
cluding T cell costimulation. 1 Naive T cells require 2 signals to 
become activated. 2 The first signal results from the interaction 
between antigen-specific T-cell receptors and an antigenic pep- 
tide presented in the context of a major histocompatibility com- 
plex class on the surface of antigen-presenting cells (APCs). The 
second set of signals is delivered via costimulatory molecules 
that are expressed on the cell surface of APCs and cytokines that 
are either produced by APCs, and/or by the activated T cells 
themselves. Costimulatory molecules belong either to the im- 
munoglobulin family members (e.g., CD28/B7) or the tumor 
necrosis factor superfamily (e.g., CD40/CD40L). 3 

It is well known that the disruption of immune tolerance con- 
tributes to the development of allergic asthma (AA) and rheu- 
matoid arthritis (RA). In this context, it is plausible that costim- 
ulatory molecules which initiate immune responses play im- 
portant roles in the pathogeneses of these diseases. 4,5 In the 



present study, we hypothesized that a common genetic interac- 
tion involving costimulatory molecules exists which renders 
susceptibility to both AA and RA. To test this hypothesis, we 
sought to identify genetic interactions that involved costimula- 
tory molecules in patients with AA using multifactor dimen- 
sionality reduction (MDR) analysis, and then to confirm the rel- 
evance of such interactions in patients with RA. 

MATERIALS AND METHODS 

Study subjects and phenotype definitions 

Patients with AA or RA were recruited at the Seoul National 
University Hospital. A diagnosis of asthma was made when a 
subject with symptoms of dyspnea or wheezing showed revers- 
ible airway obstruction in accordance with the guidelines issued 
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Table 1. Genes and SNPs evaluated in this study 



Gene(GenBanklD) 


SNP 


rs number 


Minor allele 

frpm ipnp\/ 

1 1 CLjUCI IL>y 


C7M4* (1493) 


-3180T 


rs5742909 


12.9% 




q a n 
+4aA > b 


rc 9Ql 77c 
iSZo 1 //□ 


6 1 .0 /o 




+6230G>A 


rs3087243 


17.8% 


CD28W 


IVS3+17T>C 


rs31 16496 


7.5% 


CD86[m) 


-3479T>G 


rs2715267 


24.6% 




1179V 


rs2681417 


0.3% 


CD40(%S) 


-1C>T 


rs1 883832 


30.8% 


CD40LG(%9) 


-3458A>G 


rs3092952 


6.0% 



^Cytotoxic T-lymphocyte-associated protein 4. In the control group. 



by the National Institutes of Health. 6 All asthma patients showed 
a positive skin prick test response (allergen/histamine ratio > 1.0 
plus a mean wheal size >4 mm) to 1 or more aeroallergens. RA 
was diagnosed according to the 1987 revised criteria issued by 
the American College of Rheumatology (formerly the Ameri- 
can Rheumatism Association). 7 Control subjects were recruited 
from among healthy individuals. No control subjects com- 
plained of respiratory or joint symptoms as determined during 
interview with an allergy and rheumatology specialist, using a 
supporting questionnaire. In addition, they showed negative 
results in methacholine bronchial provocation tests (PC20 > 16 
mg/mL) and in skin prick tests with common aeroallergens. All 
subjects enrolled in this study provided written informed con- 
sent, and the study protocol was approved by the Seoul Nation- 
al University Hospital Institutional review board (H- 1004-033- 
315). 

Selection of genes and SNPs 

Based on previous reviews, 4,8 5 candidate genes encoding co- 
stimulatory molecules, and 8 SNPs in those genes were selected 
(Table 1) namely, -3180T, +49A>G, and 6230OA in CTLA4 
(coding for cytotoxic T-lymphocyte-associated protein 4), IVS3+ 
17T>C in CD28 (coding for CD28), -3479T>G and I179V in 
CD86 (coding for CD86 [B7-2]), -1C>T in CD40 (coding for 
CD40), and -3458A> G in CD40LG (coding for CD40 ligand). All 
SNPs selected were functionally relevant and previous reports 
showed significant associations between them and various im- 
munologic diseases (except IVS+17T>C in CD28, see Supple- 
mental Tables 1 and 2). Scoring was performed by the high- 
throughput single base-pair extension method (SNP- IT™ assay) 
using an SNPstream25K system, which was customized to au- 
tomatically genotype DNA samples in 384-well plates and pro- 
vide a colorimetric readout (Orchid Biosciences, New lersey, 
USA). 

Statistical analysis 

Association analysis was performed for each SNP using a mi- 
nor allele dominant model {i.e., AA+AB vs BB, where A and B 
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Table 2. Summary of MDR findings 


No. of 
factor 


Combination 


CV* 
consistency 


Prediction 
error (%) 


Pvalue 


2 


-3479T>G(CTO), 
-3458A>G(CZMK) 


10/10 


27.4 


0.001 


3 


-3479T>G(CD#6), 
-3458A>G (CD40LG), 
+49A>G(C7L44) 


8/10 


32.2 


0.054 


4 


-3479T>G(C000), 
-3458A>G(CZMM), 

jTICIA S.P \PT\ A/l\ 

-]C>J(CD40) 


6/10 


36.7 


0.377 


5 


-3479T>G(CTO), 
-3458A>G(CMK), 

+49A>G(C7L44), 
-]C>T(CD40), 

-318C>T(C7L44) 


4/10 


40.3 


0.828 



"Cross validation. 



are major alleles, and minor alleles, respectively). P values and 
odds ratios were determined by logistic regression analysis con- 
trolled for age and gender. Compliance with the Hardy- Wein- 
berg equilibrium was determined using 2X2 tests. MDR analy- 
sis (version 3.0.2; Computational Genetics Laboratory, Dart- 
mouth Medical School, Hanover, NH; http://www.epistasis. 
org) was used to reveal interactive genetic effects on AA or RA, 
as described previously. 9 -3458A>G in CD40LG is located on 
the X chromosome and thus male subjects had only one allele 
(A/- or B/-). Because only 3 genotypes (A/A, A/B, and B/B) for 
one locus are allowed in MDR analysis, 10 we treated A/- as A/A 
(or B/- as B/B). Significance of SNP combinations obtained 
from MDR analysis was confirmed using a contingency table 
and a chi-square test. All statistical analyses were performed 
using the R software (version 2.15.3, www.r-project.org), and P 
values > 0.05 were regarded as significant. 

RESULTS 

A total of 200 patients with AA (age 48.4+16.7 years [mean± 
SD]) and 184 patients with RA (age 52.6 ± 12.3 years [mean ± SD]) 
participated in the study. The proportion of male was 32.0% for 
patients with AA and 24.5% for patients with RA, respectively. 
One hundred and 82 subjects were enrolled as controls (age 
47.6+13.8 years [mean+SD]) and 37.4% of them were males. 
All eight SNPs examined were in Hardy- Weinberg equilibrium. 
Logistic regression analysis, adjusted for age and gender, re- 
vealed no significant association between any individual SNP 
and AA, and between individual SNPs and RA in a minor allele 
dominant model (data not shown). We then performed MDR 
analysis in the control and AA subjects. Table 2 summarizes the 
average cross validation consistency, the average prediction er- 
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Table 3. Combined genetic effects COT-3479T>G and CD40LG -34bQA>G 
on allergic asthma 



Genotype 


Frequency(N,%) 


CD86 


CD40LG 


Control 


AA 


TT 


AA 


15(8.2%) 


53(26.4%) 


TT 


A- 


69(37.9%) 


40(19.9%) 


TT 


AG 


6(3.3%) 


17(8.5%) 


TT 


G- 


7 (3.8%) 


6(3.0%) 


TP 

lb 


A A 

AA 


1 r in oo/ \ 

lb (o.o%) 


OO /1 O 00/ \ 

zo (13.9%) 


TG 


A- 


52 (28.6%) 


28(13.9%) 


TG 


AG 


5(2.7%) 


9(4.5%) 


TG 


G- 


1 (0.5%) 


0(0%) 


TG 


GG 


3(1.6%) 


5(2.5%) 


GG 


AA 


2(1.1%) 


6(3.0%) 


GG 


A- 


4(2.2%) 


5(2.5%) 


GG 


AG 


1 (0.5%) 


0(0%) 


GG 


G- 


1 (0.5%) 


3(1.5%) 



*Bold characters represent the most frequent SNP combination in control and 
allergic asthma (AA) subjects respectively. 



ror, and the P value obtained by MDR analysis for each of the 
factors evaluated. A single 2-factor model (-3479T>G CD86 
and -3458A> G CD40LG) had a maximum cross validation con- 
sistency (10 out of 10) and a minimum prediction error (27.4%, 
P= 0.001). We re-confirmed these results using a contingency 
table and a chi-square test. We took a simple allele counting ap- 
proach to combine the results for the male and female samples 
instead of X chromosome analysis stratified by gender, which 
leads to a loss in power due to stratification. 11 We considered 
each A/ A, A/B, B/B, A/-, and B/ - as separate genotypes, and 
then combined them. Table 3 showed combined genetic effects 
of -3479T > G CD86 and -3458A > G CD40LG on AA. The frequen- 
cies of each SNP combination were significantly different be- 
tween the control and AA groups (df=12, chi-square =50.982, 
P<0.001). Subjects with a T/T genotype of -3479T>G CD86 
and an A/A genotype of -3458A>G CD40LG were found to be 
significantly more likely to develop AA than those with a T/T 
genotype of -3479T>G CD86 and an A/- genotype of -3458A>G 
CD40LG (adjusted OR, 6.09; 95% CI, 2.89-12.98; logistic regres- 
sion analysis controlled by age). 

Next, we evaluated the combined effects of -3479T>G CD86 
and -3458A> G CD40LG on RA (Table 4). Similarly, the frequen- 
cies of each SNP combination were significantly different be- 
tween the control and RA groups (df= 12, chi-square = 148.925, 
P<0.001). In addition, we found that the T/T genotype of 
-3479T>G CD86 and the A/A genotype of -3458A>G CD40LG 
were interactively related with a significant risk of developing 
RA compared to a combination of a T/T genotype of -3479T> G 
CD86 and an A/- genotype of -3458A> G CD40LG (adjusted OR, 
39.35; 95% CI 15.01-107.00; logistic regression analysis con- 



Table 4. Combined genetic effects C000-3479T>G and CD40LG -3458A>G 
on rheumatoid arthritis 



Genotype Frequency (N, %) 



CD86 
(-3479T>G) 


CD40LG 
(-3458A>G) 


Control 


RA 


TT 


AA 


15(8.2%) 


77 (41.8%) 


TT 


A- 


69(37.9%) 


9(4.9%) 


TT 


AG 


6(3.3%) 


18(9.8%) 


TT 


G- 


7(3.8%) 


3(1.6%) 


TT 


GG 


0(0%) 


3(1.6%) 


TG 


AA 


16(8.8%) 


32(17.4%) 


TG 


A- 


52(28.6%) 


11(6.0%) 


TG 


AG 


5(2.7%) 


13(7.1%) 


TG 


G- 


1 (0.5%) 


3(1.6%) 


TG 


GG 


3(1.6%) 


1 (0.5%) 


GG 


AA 


2(1.1%) 


10(5.4%) 


GG 


A- 


4(2.2%) 


1 (0.5%) 


GG 


AG 


1 (0.5%) 


3(1.6%) 


GG 


G- 


1 (0.5%) 


0(0%) 



*Bold characters represent the most frequent SNP combination in control and 
rheumatoid arthritis (RA) subjects respectively. 



trolled by age). MDR analysis in the control and RA subjects also 
revealed a significant interaction between -3479T>G CD86 and 
-3458A>G CD40LG (Table S2). 

DISCUSSION 

Immune tolerance is related to conditions such as allergy, au- 
toimmunity, tumor tolerance, and organ transplantation rejec- 
tion. A previous large population-based study revealed that Thl- 
and Th 2-mediated diseases are significantly linked, 12 which 
supports the theory that autoimmune and atopic diseases share 
risk factors that increase the propensity of the immune system 
to generate inappropriate responses to non-pathological anti- 
gens. AA and RA are examples of immune-tolerance-related 
diseases, and it is known that antigen-specific T cells play an 
important role in the development of these diseases. Because 
costimulatory signals are critical for optimal T-cell activation, 
proliferation, and differentiation, an enormous number of ex- 
perimental and pre-clinical studies have examined the possible 
use of costimulatory molecules as therapeutic targets in AA and 
RA. Interestingly, it has been reported that RNAi-mediated 
CD40-CD40LG interruption not only promotes tolerance in au- 
toimmune arthritis 13 but also significantly reduces nasal aller- 
gic symptoms. 14 Accordingly, genetic susceptibility involving 
costimulatory molecules may be a common risk factor for AA 
and RA. 

An initial MDR analysis revealed that a genetic interaction be- 
tween CD86 and CD40LG conferred susceptibility to AA and 
RA. CD40LG is located on the X chromosome. Most loci on the 
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X chromosome of females are subject to X-chromosome inacti- 
vation (only 1 allele from each pair of alleles is expressed). 15 
However, a recent report raised the possibility that some loci on 
the X chromosome of females can escape X-chromosome inac- 
tivation. 16 For example, females had higher levels of CD40LG 
than males 16 and therefore the genetic risk to a male hemizy- 
gote (A/-) may not be equivalent to the genetic risk to a female 
homozygote (A/ A), which requires alternate hypotheses for 
testing association on the X-chromosome loci. 1117 To address 
this, we used a simple allele counting approach to combine the 
results for the male and female samples, instead of X-chromo- 
some analysis stratified by the gender, which leads to a loss in 
power due to stratification. 11 It seems that the A/A genotype of 
CD40LG in a female increases the risk of AA or RA, whereas the 
A/- genotype of CD40LG in a male plays a protective role. How- 
ever, these associations become more prominent if they are 
combined with the T/T genotype of -3479T>G CD86. At pres- 
ent, we cannot fully exclude the possibility that gender-specific 
traits, other than genotypes of CD40LG, may contribute to these 
associations. 

Of the molecules involved in co -stimulatory signaling, CD40LG 
and its receptor CD40 have been shown to be one of the most 
critical. Furthermore the absence of the CD40LG-CD40 interac- 
tion results in defective T-cell-dependent immune respons- 
es. 1819 However, protective CD4 T cell responses to viral infec- 
tions found in CD40- or CD40LG- deficient mice suggest the 
existence of an alternative pathway for priming of T cells. 20,21 
The up-regulation of CD86 (B7-2) on APCs has presented a 
CD40LG-CD40-independent pathway for T-cell priming. 22 
These findings suggest that the CD40LG-CD40 pathway and 
CD20-B7 pathway are interactively involved in the immune re- 
sponse. In addition, CD86 -3479T> G was significantly associat- 
ed with atopy and an increased risk of asthma when combined 
with a promoter polymorphism (rs 1800872) in the IL10. 23 Al- 
though no direct interaction between CD86 and CD40LG has 
been reported to date, it is possible that allelic variants of the 
CD86 and CD40LG genes interact to promote the development 
of AA and RA. However, no AA patient enrolled in this study 
complained of arthritis symptoms, and similarly no RA patient 
complained of respiratory symptoms. It appears that genetic 
susceptibility involving costimulatory molecules contributes 
only to the early stages of AA or RA development, that is, the 
initiation of the immune response. Accordingly, AA- or RA-spe- 
cific immune responses following breakdown of immune toler- 
ance, can be initiated only in the presence of a genetic suscepti- 
bility or an environmental trigger. 

SNP allele frequencies are different according to ethnicity. For 
example, the minor allele frequency of CD86 -3479T>G in the 
dbSNP database shows wide variations; 30.0% in European 
populations, 26.7% in Chinese populations, 21.1% in Japanese 
population, and 70.8% in African populations. 24 In this study, 
the minor allele frequency of CD86 -3479T>G in the control 



subjects was 24.6%, which is similar to that of Chinese and Jap- 
anese populations. These data enhance the reliability of the sta- 
tistical analyses performed in this study. 

To evaluate the functional relevancies of CD86 -3479T > G SNP 
and CD40LG -3458A>G, we searched their recognition motifs 
in the 'RegulomeDB' database (http://regulome.stanford.edu/ 
index). RegulomeDB includes high- throughput experimental 
data sets from various sources, as well as computational predic- 
tions and manual annotations to identify putative regulatory 
potential, and identify functional variants. 25 Predicted binding 
motifs for -3479T > G CD86 are SRF, ELF3, TCF3, TCF7, TCF7L2, 
and ZFP105 and predicted binding motifs for -3458A> G CD40LG 
are ZFP105 and ZSCAN4, which suggests that both SNPs may 
be functionally competent. 

Before generalizing our results, 2 points should be considered. 
Firstly, the case subjects enrolled in this research probably also 
had co -morbid diseases as well as AA or RA, which could have 
confounded our results. However, the inconsistent nature of co- 
morbid diseases, and the comparison made with apparently 
healthy controls might dilute any genetic effects originating 
from the co -morbid diseases. Secondly, we could not complete- 
ly remove the risk that spurious interactions might have been 
identified, despite confirmation of our results by MDR analysis 
based on permutation. Accordingly, we suggest that our find- 
ings be confirmed by a large-scale study at a later date. 

In summary, this study shows that allelic variants in the CD86 
and CD40LG genes interact to significantly increase the likeli- 
hood of both AA and RA development. As mentioned above, 
we believe that this interaction does not entirely explain AA or 
RA development. Nevertheless, our results encourage us to ex- 
plore genetic susceptibility conferred by costimulatory mole- 
cules in the context of immune-tolerance-related diseases. 
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